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Abstract 

A quantum field theory of the liquid-glass transition in water based on the two 
band model in the harmonic potential approximation is presented by taking into 
account of the hydrogen bonding effect and the polarization effect. The sound 
and diffusion associated with intra-band density fluctuations, and the phonons and 
viscocity associated with inter-band density fluctuations are calculated. The Kauz- 
mann paradox on the Kauzmann's entropy crisis and the Vogel-Tamman-Fulcher 
(VTF) law on the relaxation times and the transport coefficients are elucidated 
from the sound instability at a reciprocal particle distance corresponding a hydrogen 
bond length and at the sound instability temperature very close to the Kauzmann 
temperature. The gap of specific heat at the glass transition temperature and the 
boson peaks are also presented. 
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1 Introduction 



Water plays an essential role in a medium in which chemical and biological molecules 
operate and it influences their functions. Its properties are far from being fully under- 
stood. Experimental studies of the glassy behaviours of bulk water are extremely difficult, 
while hydration water particularly on the surface of protains can provides an interesting 
alternative to study the slow dynamics of water, because the binding of water molecules to 
the surface can prevent crystallization. Experimental studies in hydration water certainly 
show the universal features in the liquid-glass transition [TJ |2l El H] • We intend to develop 
theory of the liquid-glass in hydration water. Here at the first stage we develop the glass 
transition in bulk water. 

We have established quantum field theory of the liquid-glass transition in a one- 
component liquid in the framework of the two band model to elucidate the Kauzmann 
paradox, the VTF law, the gap of specific heat at the glass transition temperature and 
boson peaks E] , which appear universally in the liquid-glass transition. We have also 
proposed a quantum field theory of phonons in a two-component liquid [5J |7J |S], and 
also proposed a quantum field theory of phonons and melting in compounds taking into 
account of the electric polarization in the rigid ion approximation, the Ewald method 
[3 EH E] • Here we extend the quantum field theory of the liquid-glass transition in a 
one-component liquid to that in water considering the quantum field theory of phonons 
in a multi-component liquid and the electric polarization in the rigid ion approximation 
[3 HII1HI]- The essential points of our theory are as follows 0E]: 

(1) The two band model. A particle is temporarily in a randomly distributed harmonic 
potential made up by the surrounding particles, making up and down transitions between 
the ground state and the first excited states, and then hops to a surrounding vacancy. This 
picture constitutes the model Hamiltonian which consists of the random eigenfrequencies 
and random hopping magnitudes. The configurational averaged model Hamiltonian yields 
an unperturbed Hamiltonian. Due to hopping, instead of individual states of the local 
potentials, the system have energy bands having a width determined by the hopping 
magnitude. 

(2) Intra-band and inter-band density fluctuations. Density fluctuations consist of 
intra-band density fluctuations related to the hopping to a neighbouring vacancy and 
inter-band density fluctuations associated to the up and down transitions between two 
bands. Intra-band and inter-band density fluctuations are associated with intra-band and 
inter-band elementary excitations, respectively. The collective excitations of intra-band 
and inter-band elementary excitations are sound and phonons, respectively. Sound is a 
density fluctuation wave. In the classical limit the particle energy dispersion associated to 
intra-band elementary excitations tends to a free particle one so that the sound reduces 
to a density fluctuation wave .12,. Phonons are elastic waves. In the classical limit 
coupling strengths between particles deviated from their assigned positions constitute 
spring constants so that the phonons reduce to elastic waves [T2*] . 

(3) Dissipation and relaxation times. Random scattering processes come from the 
random potentials which is described by subtracting the unperturbed Hamiltonian from 
the model Hamiltonian. Dissipative process comes from the simultaneous scattering pro- 
cesses of two particles by the same random potential. The sum of the configurationally 
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averaged simultaneous scattering processes over all sites yields the correlation functions 
of random harmonic frequencies and random hopping magnitudes, which leads to the 
relaxation times; the former and the latter correspond to the f3 and a-relaxations, respec- 
tively. When the relaxation times of sound and phonons are shorter than their respective 
periods, diffusivity and viscosity appear, respectively. 

(4) The Kauzmann entropy and the VTF law. Intra-band density fluctuation entropy 
consists of sound entropy, fluctuation entropy due to intra-band elementary excitations 
and dissipative entropy due to diffusion. In an equilibrium state the dissipative entropy 
compensates the fluctuation entropy with a negative value and the system reaches a 
local equilibrium. But since sound depends on temperature, a mixing between the sound 
and fluctuation entropies occurs. The mixing entropy yields the Kauzmann entropy: a 
Curie law with a minus sign, where the critical temperature corresponds to the sound 
instability temperature which indicates fragility. The hopping amplitude relates to the 
multi-dimensional configuration space of the Kauzmann entropy; the hopping amplitude 
is proportional to the exponent of the Kauzmann entropy per particle. This hopping 
amplitude explains the VTF law: the exponent of the Curie law with a minus sign per 
particle. The hopping amplitude is the origin of the a-relaxation. 

(5) The glass transition temperature and the magnitude of the randomness of harmonic 
frequencies. Inter-band density fluctuation entropy also consists of phonon entropy, fluc- 
tuation entropy due to inter-band elementary excitations and dissipative entropy due to 
viscosity. In an equilibrium state the dissipative entropy compensates the fluctuation 
entropy with a negative value and the system reaches a local equilibrium. Since phonons 
do not depend on temperature, mixing does not occur. The magnitude of randomness of 
harmonic frequencies relates to the mult i- dimensional configurational space of the fluctu- 
ation entropy; the magnitude of randomness of harmonic frequencies is proportional to 
the exponent of the fluctuation entropy per particle. The fluctuation entropy prevents 
the Kauzmann entropy crisis at the glass transition temperature. The magnitude of the 
randomness of harmonic frequencies is the origin of the /3-relaxation. 

In applying the two band model to the liquid-glass transition in water G2] , we 

should clarify the structure of water. The characteristic features of the structure of water 
are as follows (i) Hydrogen bonds play an essential role in clustering water molecules 
comparing with the common molecules. A water molecule consists of an oxygen and two 
covalent bonded hydrogens so that the water molecule has the two protons and the lone 
pairs. A proton of a water molecule is hydrogen bonded to a lone pair of the nearest 
neighbour water molecule, (ii) Water has the proper geometric structure. Water is made 
up of a strongly hydrogen bonded network and locally very similar to ice. But these 
hydrogen bonds have very short lifetimes, (iii) A water molecule is electrically polarized, 
which originates from the hydrogen covalent bonds. 

In order to include the above characteristic features, first we determine the interac- 
tion potentials V a p with the corresponding pair distribution functions g a p, which predict 
the structure of water, where a, (3 mean an oxygen atom O or a hydrogen atom H. (i) 
The interaction potential between O-H, Voh consists of hydrogen covalent bonding, Vq^ 
and hydrogen bonding, Vq^. A strong hydrogen bonded network is constructed by the 
couplings of VoH b V^o an d Kdh^ho 3 ) which comes from the fact that the configurational 
average of the term VqhVho over the position of the hydrogen atoms yields the terms 
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VJh Vho and Wo. not the terms v oh v ho and VoMo> because the former terms 
are favorable. The hydrogen bonding plays an essential role in clustering water. The 
first peak of gou(R) corresponds to the hydrogen bond length [T^j. Since hydrogen bond- 
ing plays an essential role in water, the sound instability occurs at a reciprocal particle 
distance corresponding to a hydrogen bond length, (ii) In the interaction potential Vhh 
we take into account of the interaction between hydrogens in different water molecules. 
Hopping between protons and lone pairs reflects Vhh and the hopping yields very short 
lifetimes of hydrogen bondings, (iii) There exists the electric polarization between O-H. 
The electric polarization relates to the p-state of particles. An excited state of a parti- 
cle corresponds to a p-state, while a ground state to an s-state. Thus we need not take 
into account of the polarization in the interaction potentials between intra-band density 
fluctuations because the intra-band density fluctuations annihilate and create the same 
states. However we should take into account of the polarization between inter-band den- 
sity fluctuations because the inter-band density fluctuations annihilate an s-state and 
create a p-state, and vice versa. We take into account of the electric polarization in the 
rigid ion approximation, the Ewald method |§1 ITUl ITT]. 

We organize as follows: in section 2, we briefly survey the formulation. We obtain the 
model Hamiltonian with random harmonic frequencies and random hopping magnitudes. 
Taking the configurational average of the model Hamiltonian, We obtain the unpurterbed 
Hamiltonian. We construct the interaction Hamiltonian for intra-band and inter-band 
density fluctuations, which leads to the sound and the phonon modes. In section 3, we 
introduce the interaction Hamiltonian for the random scattering processes, which leads 
to the dissipative process. In section 4, the correlation fuctions of intra-band density fluc- 
tuations with the vertex corrected bubble diagrams of intra-band elementary excitations 
yield the sound and diffusivity. The sound instability occurs at a temperature T and at a 
reciprocal particle distance corresponding to a hydrogen bond length. In section 5, taking 
into account of the electric polarization due to hyrogen covalent bond, we calculate the 
correlation functions of inter-band density fluctuations with the vertex corrected bubble 
diagrams, which yield phonons and viscosity modified by the polarization. We also show 
the boson peaks. In section 6, we obtain the Kauzmann entropy from the mixing of 
the sound and the intra-band fluctuation entropies in the sound instability regime. We 
show that a unique hopping amplitude corresponding to the Kauzmann entropy governs 
the velocities of individual particles, the sound and the hopping amplitudes of individual 
particles, while the magnitude of the randomness of the hamonic frequencies of particles 
depends on the individual particles. The glass transition occurs at the crossover temper- 
ature of the Kauzmann entropy and the intrer-band fluctuation entropy. We derive the 
specific heat and show the VTF law for the relaxation times for sound and phonons, and 
the transport coefficients. Section 7 is devoted by some concluding remarks. 

2 Formulation 



Here we briefly survey the formulation. We start with the following Hamiltonian: 
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where Tie° a (— iV) is the energy operator of a free a-particle, n a (x) = ip a (x)ip a (x) . We 
write the operator ip a (x) by using the localized operators, which are rewritten by the 
extended operaters. 

1p a {x) ^ ^ {x Ram)b a mu- b ama , ^ ^ 6 ^ 0-afj.p- (2-2) 

The Heisenberg equation for 6 QmAt is given by 
d f 

-h—b amfM = ^2 / d 3 xw amu (x - R am ) [he° a (-iV) + ^^(x)] w anu {x - R an )b anv , 

nu 

(2.3) 

= J2 I d 3 yV a p(x - y)n^y), (2.4) 

P J {Ran} 

where {R an } means summing the surrounding sites of the particle at R an in np(y). The 
term ^n an means that a potential is determined by the surrounding particles of the site 
Ran- Thus we obtain the Schrodinger equation: 

(fte°(-iV) + Qrv^x)} w ani/ (x - R an ) = huj anv w anu (x - R an ), (2.5) 

where u anfM and w anu are the eigenfrequency and the eigenfunction of the /ith state at the 
site R an in a simmilar manner to a one-component liquid in the harmonic approximation. 
We consider the following relations: 



W n al{%i) — — 2CrtaVjU> raQ ,o( a; i); Wnao( x i) ~ ^Cna^ iW n ao{Xi) ] C 



2M a U na ' 



(2.6) 

where C, na is the mean width of the zero point motion, M a is the mass of the a-particle 
and u na is the harmonic frequency; Tiuj na o = \fiuj na and Tiuj na i = \hoj na - At this stage, we 
can take the following model Hamiltonian in terms of localized operators: 

H = ^ ] huiamfjibamubama + ^ ] ^ J am^anvb\ mj ]3 a nv ■, (2-7) 
amu amfinv 

where the hopping matrix is given by 

^"Jamfianu J d XW a mui,X R a rri) ^Ra n (^Xjw anu {x R an )- (2-8) 

Note that the interaction Hamiltonian is nonlinear through the term np{y) in $h q „ and 
the potential ^R an permits a particle to hop only to a vacancy. 

Now we start with the new unperturbed Hamiltonian which is the configurational 
average of the model Hamiltonian (2.7): 

H =< H > c = ^2 ^a^L^am/i + hJatiiRam ~ Ran)b amfM b an u, (2.9) 



afim aumn 



where J afl is the hopping magnitude. < > c means the configurational average. The 
configurationally averaged values do not include the indices of the positions. Here for 
simplicity we limit ourselves to the hopping between like atoms and the hopping to those 
between the same levels. Using Eq.(2.2), we obtain 



' ' -'•a 

a^p m 

(2.10) 

From Eq.(2.2), we obtain 

n a ( X ) = J2 e ~ i ' 1 ' B<Xm ^2^ X ~ R <*™)pl dm -^{aVAix ~ Ram)pl iq + '••], 

"a 

mq /i i 

(2.11) 

where 

Padfiq — y ] a au V a aiJ,p+q, Paiq = ( a ai v a aOp+q + a aOp a <MP+<j) • (2-12) 

p p 

Pad/xq an d Pmq are intra-band and inter-band density fluctuation operators, respectively. 
Substituting Eq.(2.11) into the interaction Hamiltonian in Eq.(2.1) we obtain the inter- 
action Hamiltonian: 

a/3[ii/q V a P aflijq V a P 

Kfi,(q) = v^W. V&jiQ) = VnJWmm, (2.14) 



x y d 3 xd 3 yw a ^(x - R am )w afl ,(x - R am )V a p(x - y)w Pu (y - Rpn)wp v >{y - Rpn), 

(2.15) 

where the prime on X means if/ excludes loop diagrams, which are taken into accout in 
the potential ^R an . 

3 Dissipative process 

In randomly distributed particles, we have a model Hamiltonian, Eq. (2.7). Tak- 
ing the configurationally averaged Hamiltonian of Eq.(2.7), we obtain the unperturbed 
Hamiltonian, Eq.(2.9) and (2.10), which constitutes the two bands. The intra-band and 
inter-band elementary excitations are excited around the two bands and they are associ- 
ated with intra-band and inter-band density fluctuations, Eq.(2.12), which constitute the 
model Hamiltonian (2.13). We employ the model Hamiltonian (2.13) to describe the dy- 
namical process. The collective excitations for the intra-band and inter-band elementary 
excitaions are sound and phonons. The life times of sound and phonons appears when the 
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dispersion curves merge into the continuum of the intra-band and inter-band elementary 
excitations, respectively. 

Here we investigate the dissipative process from the random scatterings due to the 
random harmonic frequencies and the random hopping magnitudes. The original model 
Hamiltonian, Eq.(2.7) reflects the randomly distributed particles, while the unperturbed 
Hamiltonian, Eqs.(2.9) and (2.10) reflect the free particle picture. The random scattering 
process is represented by the interaction Hamiltonian Hf. 

Hi ^ ^ h{u am ^ ^a^i)^ am ^b am ^ ^ ^ ^{Jam^ian^i Ja^Rctm Ran)}b\ im ^b an f JL . (3.1) 
amfi amnfi 

Note that since the term J am ^an^ involves operators n(y), Hi is essentially nonlinear. 
But at low temperatures an atom stays longer at the same site. Thus we make the 
approximation that the term J am ^ Qn ^ is a random c-number. Using Eq.(2.2), we rewrite 
Eq.(3.1) as 



m 

+ \^fllJ - J (R -R Up-iP'-Ram+ip'-iiar. 



a Lp a ^ P '- (3-2) 



Dissipative process is constructed by elementary scattering processes due to simulta- 
neously scattering processes of two particles in /i and uth bands by the same random 
potentials. Elementary scattering processes come from the scattering processes due to 
random harmonic frequencies and random hopping magnitudes, Eqs.(3.1) and (3.2). The 
sum of configurationally averaged elementary scattering process over all sites due to ran- 
dom harmonic frequncies yields a correlation functions of random harmonic frequencies, 
U^ u and that due to random hopping magnitudes yields a correlation function of random 
hopping magnitudes, Uj^, which are given by 

^w' 1 = "at2 ^ ] (^ &a[an&avn >c ~~ ^ olv) i (3-3) 
a m 

Uj^ (<?) = Jy- J d Rg(R)e 9 {< Jam^nfM-Jamunu >c ~ ^ J an{R) J cw(-R)}| H=R ctm -H an ■ 

(3.4) 

Eqs.(3.3, 3.4) yield the relaxation times; the /3 and a-relaxation, respectively. The elemen- 
tary scattering processes, U" correspond to that in the case of two electrons scattered 
simultaneously by the same random impurity in electric conductivity. The statistical av- 
erage on a random system constrains the dynamical processes, which leads to dissipation. 



4 Sound and diffusivity 

We first investigate the sound Green's function. In the low temperature regime, we 
confine ourselves to the lower band. 

V a @ iu„ 
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where we should not confuse the temperature j3 = 1/fcsT and the suffices j3\ the temper- 
ature (3 appears as an independent parameter with the inverse energy dimension, while 
the suffices (5 stand for O or H. In the random phase approximation, we obtain 



F a p{q) = P a {q)5 a p + P a (q)V ai (q)F^, (4.2) 



?o(go + £)-«) 2: 



PM = 7 ,\T , n „ , (4-3) 



(4.4) 



\ l\ a l\g , 
V " M mf=n 

= ^ J d 3 Rg a(3 (R)e^ R V aP (R) 

V a p{R am - Rf3n) = J d 3 xd 3 yw 2 am (x - R am )V a p(x - y)w 2 pn (y - Rp n ), (4.5) 

^- = jN a (US°° + Uf ) (4.6) 



(<) 2 = ^<z 2 , ^ P 4E%. ( 4 - 7 ) 



N ^ dp 



p 



(4.8) 



where we have analytically continued \v n — > g i n Eq.(4.3). r a o and the relaxation 

time and the mean velocity of the a-particle. Eq.(4.3) is presented in Hereafter, we 
omit the suffices d and in V^ 0/30 . 

Appling Eq.(4.2) to water, we obtain the dynamical equation: 

1 — PoVoo PoVoh \ ( Poo Pqb\ _ ( Pq 
P n V U o l-PnVun)\Fuo F nu J~ \0 P H 

with the solution: 

Foo F OB \ 1 f(^~ -Ph^hh)-Po PoVohPh \ q \ 

F KO F nn ) det | |V ^h^ho^Ph (1 - PoV Q o)Pn) ' { ' 

det | | = (1 - PoVbo)(l - PhVhr) - PoVqhPhVho- (4.10) 

Here in order to construct the hydrogen bonding network, we should consider that Von(R) 
in Eq.(4.4) consists of hydrogen covalent bonding, Vq^(R) and hydrogen bonding, Vq^(R), 
the pair distribution functions of which are g^(R) and g^(R), respectively. Correspond- 
ing to this fact, we consider that Voh(?) consists of hydrogen covalent bonding VQ^(q) 
and hydrogen bonding V^(q). Thus we put 

VbH(g) = V^(q), V^(q). (4.11) 

But in deriving Eq.(4.8), we have used Eq.(4.4) under the nearest neighbour approxima- 
tion so that the terms in Eqs.(4.9, 10) involves unphysical cross terms such as Vqh Vjjq 
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and Vqh^ho- When we calculate Eq.(4.1), we should configurationally average all pos- 
sible Feynmann diagrams over the position of particles R am - If we configurationally 
average the term < VohVho > over Rum considering the pair distribution functions 
9oh(R) an< ^ PokC-^)' the mos t dominant terms in < Voh^ho >c come from the terms 
VoH C ( a oH)^Ho( a Ho) anc ^ K3H( a OH)^HO c ( a Ho)> where aQ^' hb are the mean distances of hy- 
drogen covalent and hydrogen bond, respectively. Since in our calculation we concentrate 
ourselves to the hydrogen bonding network, we should take into account of both the terms 
Vq^V^q and Vq^Veo as P arr s. Thus when the term Voh^ho appears, VorVuo should be 
replaced by V^V^ + V&V£%. 

We can obtain the modes by solving the secular equation of Eq. (4.8): 

{1 - P (q)V 00 (q)}{l ~ ^h(?)^hh(?)} - Po(q)Pu(q)V OK (q)Vuo(q) = 0. (4.12) 

Before we solve Eq.(4.12), we investigate the secular equation for the sound in a system 
consist of only a-type of particles: 

1 + P a (q)v aa (q) = V (4-13) 

We investigate Eq.(4.13) in the two limiting cases: 
(i) qoTao > 1 

Eq.(4.13) leads to the sound frequency for a-particles: 



q 2 = = (1 +W*M)K q ? = VTfv = *Lt, (4-14) 



«) 2 
■ ' S a (q) 



where S a (q) and v as are the static structure factor and the sound velocity of the a-particle: 

S a (q) = (1 + PV aa (q)r\ (4.15) 

- = t5<o)' (4 - 16) 

Next we investigate the sound instability. In the first aproximation, if we consider a 
first peak for a pair distribution function, we can put 

9aa(R) = Pa5(R - a QQ ), (4.17) 

where a aa is the mean distance between the a- particles. Then we obtain 

T T / \ A 2 / , S1H d'OLOL^ / \ 

Vaa{q) = ^a aa p a V aa {a aa ) , (4.18) 



aa — 



where V aa has the minimum negative value at a reciprocal particle distance K t 
37i/2a aa . Since 

lim P a (0) = -p, (4.19) 
qo-*o 

and V aa has the minimum negative value at a reciprocal particle distance K aa , the sound 
instability in the a-liquid occurs at the temperature T a0 and the reciprocal particle dis- 
tance K aa in Eq.(4.15): 

I + P a0 V aa (K aa ) = 0. (4.20) 
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(ii) q T a0 < 1 

Eq.(4.13) leads to the diffusion mode 

go = ir a0 vl s q 2 , (4.21) 
where the diffusion coefficient is given by 

D a = r a0 v 2 as . (4.22) 
Next we investigate the secular equation Eq.(4.12) for water in the two limiting cases: 

(I) qoTaO > 1 

We obtain the sound modes: 
<?o = "Ho) 

= v 2 ± q 2 , (4.23) 
where v s ± are the sound velocities: 

\j \{vL + <±\/ « - <) 2 + ^ 2 v 2 o P < p Vou(q)Vno(7}- (4-24) 



If the term VohVho dominates, cj s _ and v s - disappear. Next we investigate the sound 
instability. In the first aproximation, if we consider a peak for the pair distribution 
function, we can put 

g a p(R) = p P 5{R - a a/3 ), (4.25) 
where a Q/ g is the mean distance between the a and (3 particles. Then we obtain 

V aP {q) = 4nal^V a pM S -^^, (4.26) 

a a /3q 

where V a p has the minimum negative value at a reciprocal particle distance K af3 = 
37i/2a a/3 . 

In water the interaction potential Vqh and the pair distribution function g Q \ i n ^oh are 
the most important. The peak of g b n corrsponds to the hydrogen bond, which constitute 
the cluster of hydrogen bonded molecules. If we denote the mean distance corresponding 
to the peak by a h K , the reciprocal particle distance is given by The sound instability 

occurs when the secular equation Eq.(4.12) with Eq.(4.19) becomes zero at a temperature 
T and at a wave vector K = K^ n : 

{1 + (3 V oo {k)}{l + PoVhh(K)} - $V OH (K)V HO (K) = 0. (4.27) 

T is given by 

= 7^-{Voo(K) + Vuh(K) + J (Voo(K) - V HU (K)) 2 + 4V h (K)V K o(K)}. (4.28) 
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(II) q T a0 < 1 

We obtain the diffusion modes: 



(4.29) 

where D ± are diffusion coefficients: 

D± = \{ r OQV 2 Qs + Tho^Hs ± \j ( r O0^ O s - r H0^Hs) 2 + ^^OO^Op^O^Hp^bH (?) V H o(<?) } 

= ^{Do + ± y (Do - Ai) 2 + ^ 2 r ,vl v r m vl p Von{q)V w {q)}. (4.30) 
If the term VohVho dominates, D~ disappears. 

5 Phonons, boson peaks and viscosity 

Next we investigate the phonon Green's functions, which are defined by 
Devlin - r 2 ) = - L— < T^MpI jq (r 2 ) > c = J- J2 e-^-^D aipj (q), 



(5.1) 

In the random phase approximation, we obtain the phonon Green's functions and the gap 
equation [1]: 

DaipM) = Qa{q)SafiSij + Q a (q) ^ < V^q) > c D jW (q), (5.2) 

7/ 

= Q«(0) < V^(0) > c -^=, (5.3) 

where 

C^(<?) =< 4== E ^^-^^-(fl™ - Hn) > c , (5.4) 

Ki/3j(i2am - Rpn) = J d 3 xd 3 yX7 iW 2 am0 (x - R am )V a/3 (x - y) VjWp (y - Rp„). (5.5) 

QM = - ( J^'X — 2 > (5-6) 

(g + i/2r QM ) i - Wa 



where r a u is called the Maxwell relaxation time and 

11111 



2r aM 2t q0 2r Ql 2r° 2r° x 
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(5.7) 



T a o,i are the relaxation times of the particles in the lower or upper bands, respectively. 
r° 1 comes from the simultaneously scattering processes of particles between the lower 
and upper bands. This term corresponds to the cosine-function term with the scattering 
angle of an electron in the electric conductivity. 

At low temperatures, since A a = 7^, the gap equation (5.3) leads to 

1 = CQ a (0) J2 < V ^i( R ?n) >c ■ (5.8) 

Substituting ( a =< ( an > c in Eqs.(2.6) and (5.6) into Eq.(5.8) under u a r aM <C 1, we 
obtain 

M a u 2 a = ~J2< V ™/3i( R Pn) >c • (5-9) 

Pn 

Eq.(5.9) means that every a-particle feels the same statistical averaged potential. From 
Eq.(5.4, 5) the matrix element of the dynamical equation (5.2) is written as 

M« - QM < V^{q) > c = -CQ a {q)M a [q {q + — )5 af3 5 ij - M ai ^{q)}, (5.10) 

TaM 

where M ai pj is the dynamical matrix for phonons. At low temperatures the dynamical 
matrix is given by 



MaipAq) = ^ E / d 3 Rg ay (R){5 Py ^ ^e^ R - M«> W*)- (5.11 



Here, as discussed on the term Voh(<z) in Eq.(4.11) in taking into account of the 
hydrogen bonding network, we should consider that the term MoiHj(<?) consists of the 
term MQ^{q) corresponding to the hydrogen covalent bonding and the term Mo^^q) 
to the hydrogen bonding. Hereafter, in the same manner as the sound, if the cross term 
M mj(q)M ]ii0j (q) appears, M omj (q)M m0j (q) should be replaced by M^ j (q)M^ i0j (q)+ 

Now we investigate Eq.(5.5). We separate the potential V into the short range one, 
V N and the Coulombic one, 

yC. 

Vai/3j(Ram ~ Rf3n) = ^aiPj{Rara ~ Rpn) + Vaif3j(Ram — Rpn)- (5.12) 

When we take into account of the polarization of the hydogen covalent bonds, we make the 
rigid ion approximation, the Ewald's method PIEIIE]- In the rigid ion approximation; 
u^o(-R) — > S(R), we obtain 



v ^ - -oiarM' (5 ' 13) 



where we have used the SI units, eo is the dielectric constant in vacuum and e a is the 
charge of an a-particle; eo = — 2en. We define the Fourier transformed potential as 

v ai pj(q) =< E ^ q ' (Ram ~ R " n) y^ 3 {R am - Rpn) > C = v a %j(q) + $SpM (5.u) 

n 
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Next we investigate the Coulomb potential: 
which satisfy 

<P%,(R) = (2^)3 / A^V 7 *?)- (5.16) 

In the Ewald method, we define a potential, ip^^R), which is the potential energy of a 
point charge e a at a distance i? from a three dimentional Gaussian distribution Pp(R) 
with the total charge e^: 

P G ,{R) = e,(-) 3/2 e"^ 2 , = epe-W, (5.17) 

7T 

which satisfy 

P?(R) = J^yJ Ap?*?)^* (5.18) 

Thus we obtain 

The Poisson equation: 

VVS»(i2) = (a^s / A(-9 2 ^S»(?))^ R = - (5-20) 



and Eq.(5.17) lead to 



= ^e"^. (5.21) 



Now we introduce 

= V%( R ) + - = <P%>(R) + '■Pap(R) (5-22) 

Using Eqs.(5.15), (5.19) and f™ d 3 R'p%(R') = e p , we obtain 



/•/PR 

Here we define 



(5 - 23) 



^VSbM+V^M. (5.24) 
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2 

Ft^Ram — Rr„ c 



First we investigate the term K' 7 , which is written as 

- - < z< " - / ***-vv)>i 

~ (27r) 3 e y g' 2 

=< X) / + ^|{e- 9 ' 2 / 4P - i}] e i(«'-«)-(^-^n) >c 

- f d 3 Rg a ,(R) [ d*q'[^ + ^{e~^ - l}]e^> R 



(2Tr) 3 e v a 

^Pim + m {e -<?/±P _ (5.25) 



eov a ' q 2 q 2 



where we have made the approximation; g a g(R) = pp and v a is a primitive unit cell. The 
term V H is given by 

^ (27T)3 eo R J^ R 1**--^. >c 

(27r) a e w a 7 

S / a ff^ P 3 / 2 / d 3 RHij(VPR)e~ [qR , (5.26) 
(27r) J e u a i 



<9 2 2 1 r°° , * 



^•(x) = ^— r -=- / rfse" 5 . (5.27) 

OXiOXj y/TC X J x 

It should be noted the function Hij{x) is singular at x = 0. When we treat the crystalline 
state [TQl [TT] , Hij(x) is replaced by if?- = d f. dx . ^^ Jq dse~ s2 , because we should exclude 
the singularity, since the dynamical matrix includes R m = R n . But in our case, taking 
into account of the spherical symmetry ^ Ha = -j^z~ x and the property of the argument 
of integration d 3 R = 4irR 2 dR yields the integration finite. Thus Eqs.(5.25) and (5.26) 
lead to 

vSfub) = - QiM)l (5-28) 

QiM = -^-[e-i - 1] + / d*RH l0 {^fPR)e-^ R , (5.29) 

toV a 47re f a J 

where the first term and the second term in Eq.(5.28) relate to the macroscopic electric 
field and the Lorentz field, respectively. In the spherical symmetry, we obtain 

Qij(q) = 5ijQ(q) = Sij-^—. (5.30) 
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There is the electrical charge neutrality: 



(5.31) 



Next we investigate the short range potential V^. In order to do so, we introduce the 
couplings like the chemical bondings: 

V ffa p(R) = - J d 3 xd 3 yV X3 w 2 a0 (x)V a/3 (x - y)V m w%{y - Re 3 ), (5.32) 

V vaP {R) = - J d 3 xd 3 yV Xl w 2 a0 (x)V a/3 (x - y)V yi w%(y - Re 3 ). (5.33) 

Here we take q = qe z and £ = qR. Since we are concerned with the spherical symmetry 
distribution of particles, we obtain 



M axax (q) = M ayay (q) = 



d 3 Rg aa (R) 



cos £ : sin £ - 

T 3 



V aaa (R){-^ + ^} 



+ K aa (R){— + + e aPa — 

V aa7 (R) 2V Kaj (R) 



+ 



}■ 



(5.34) 



M azaz (q) 



d 3 Rg aa {R) 



Va aa \tt){—£— H ^ ^5— j 



+ Kraal-") { 1 — - } - e a p a 



Vaay(R) 2\4 Q7 (i?) 



+ 



(5.35) 



M axf3x (q) = M ayf3y (q) = -^-^J^ J d 3 Rg a(3 (R) 



M a ( a y l\/3 

, cosf sinf. Tr /T _.sin£ cosf sin£. 2 
x <^ V aaf3 (R)( — ^ + — ^) + V;^(i?)(— i + — ^ - — ^) + e^e^ 



3e VaJ ' 
(5.36) 



2cos£ 2sin£ 2 
+ {R){ 75 1 -pr— ) ~ e a epPp 



3e ^aJ ' 



(5.37) 
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where the other terms are zero. Since we have taken q = qe%, the z-component of the 
dynamical matrix corresponds to the longitudinal mode and x and y-components to the 
transverse modes. It should be noted that the dynamical matrices have the effects of the 
electric polarization in the last term in the curly brackets. But because of the charge 
neutrality (5.31), the last terms in Eqs.(5.34) and (5.35) do not have the effect of the 
polarization. 

The secular equation of Eq.(5.2) leads to 

{<?o(<?o + —) ~ M 0t0 i{q)}{qo{qo + —) - M mm (q)} - M 0i m{q) M m0i (q) = 0, (5.38) 
Tom t~hm 

where i = z corresponds to a longitudinal mode and i = x, y to two transverse modes. 
Now we investigate the secular equation (5.38) in the two limiting cases: 

(I) qoTaM > 1: 

We obtain phonon modes 



ql = 

l 

" 2 



M 0i0t (q) + M mm (q) ± y/{MoiOi(q) ~ M mm (q)} 2 + 4M 0im (q) M m0l (q)] . 



(5.39) 



Note that the term M a i a i consists of the ath individual phonon frequency and the fre- 
quency shift due to the coupling of particles with different masses, whereas M ai ^ plays a 
role in the mixing between the individual phonon frequencies. 

Next we investigate the phonon frequencies in the two limiting cases: First we intro- 
duce 

vSS = ±-J d'Rg^R^C, vH = d*R 9aa (R) 2 -^pe, 

y ,(n) = J_ [ j3 T?n ( m VtraP tn „/(«) = J_ I /?„ .[ mlllHl 



era 



M, 



l - J d 3 Rg a ,(R)^fC, = ±-J d*Rg a ,(R) 2 -^e, (5.40) 

M a 6aPa 3e v a ' 



1 2 1 . 

c « = ^r e «^^7TT' ( 5 - 41 ) 



(i) qR = £ « 1 
If we define as 



A _ A _ A _ AO) , AO) _ r _ A _ r f> = B — R _V*a V*a 

■"-ax — s^-ax sT-ay — v aa ~ v na ^a — "-a ^ai J -'at — u ai ■ LJ ay — -^g 

(2) ..(2) 



(2) ? .(2) 
era Una 

io" ~5~' 



A , = A = v m + ?/ (0) + 2c = A +2c B , = B = — + — 

■rt-al — S±az — U aa T V na i" ZC a — /l a -|- ZC a , -D a / — -D a2 ; — 10 ' 

? /(2) ? /(2) o,./(2) /(2) 

£} a t = rS ax — £f ay = — — H — , h a i = 13 az = — — h ~~T7r~ 1 4 ^J 

1U o 1U 1U 
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we obtain 



M axax M a y a y A a t + B a i A a -\- B t 
M azaz = A a i + B a i = A a + B a i + 2c Q , 



at u a 



M a y a y = ir-\ — A at + B' at } = ^-a —^{—A a + B' a t + C a } , 



(5.43) 



In the long wavelength limit, we obtain the phonon frequencies: 



ql = u 2 t± {q) 



A + A u - co - c H + A^^iBot - B' m ) + Ao lX~-co-c^ ~ B 'ot), 



^ +l°_ 2 c c °_ CH (- B Ht + B' m ) + x^fc^3^( 5 ot + B' ot ). 



A +A H -co-cn 



(5.44) 



ql = wi±(q) 

Ao +A n + 2c + 2c H + - B' m ) + Ao+ 2 ^ + £% H (^ - 5'oi), 

(5.45) 

cjA y + are the optical modes and u)\ q - the acoustic modes. The velocities of the acoustic 
phonons are given by 

2Aq-2c 2A h -2c h ~ , ~ ,1 ( . 

v t = \ —(Bm + B' m ) + — {B 0t + B'ot)- (5, 46) 

A + A n - c - c H A + A H - c - c H g 

2Aq + 4cq ~ ~ : ' 2A H + 4c H ~ , ~ ,1 

V ' = VAo + A H + 2co + 2cH (i?m + ^ ffl) + Ao + A h + 2co + 2ch ( ^ + B °^ 

(5,47) 

(ii) gi* » 1: 

In the short wavelength limit, we obtain 



M _ yr - 7 ,(0) , (0) , /(o) , /(o) _ „ _ 

J "ct:ra:r lvl ayay v aa ~ v -KCt ~ v aa ~ v na ^a — '-'a ^ai 

M _ .,(0) , (o) , AO) , AO) , 2 _ o , 2 



M ax p x — M ayf3y — -jrrC 



M azf3z = ^J^(-2c a ), (5.48) 
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Then we obtain 

<?o = u t ±{q) = ^{S + S R - c - c H ± a/ (So - S u - c + c H ) 2 + 4c c H }, (5.49) 

g 2 = w,±(qf) = ^{S + S H + 2c G + 2c H ± a/^o - S H + 2c G - 2c H ) 2 + 16c c H }, (5.50) 

These modes correspond to the boson peaks, where the density of states has peaks. In 
nonpolarized liquids, because of c a = the corresponding logitudinal and tranverse fre- 
quencies and the boson peaks are of the same. 

Next we investigate the secular equation (5.38) in the low frequency regime: 

(II) q r aM < 1 

Eq.(5.38) leads to 

(Qo + ir aM M aiai )(q + iTpM Wi ) - \T aM M ai piiT m Mp iai = 0. (5.51) 
Thus we obtain the dissipative modes due to viscosity 

qo = -[ToMMoiOi + thmMekik - a/ {j 0M M 0i0i - r HM ^Hiffi) 2 + 4r M^fomiTHM^ffiOi], 

-T-HM(Bffi + 5 ffi ) + i j ; t— ToM{B 0i + B oO, (5.52) 



where the modes corresponding to the optical modes disappear because of the large imag- 
inary value. We should not confuse the imaginary number i and the suffices i; the i in 
front of the terms means the imaginary number, but i in the suffices means the component 
of the coordinates. Thus viscosities are given by 



n — / 2tqm(^4o-co) ~ 7~r> T D/ 'j i 2t H m(^H~c h ) ~ J~r> T Ul \ 

Vt-y roM(Ao-Co)+THM(AH-CH) rHM ^ Ht + ^ HtJ + tom (Ao- C o)+THm(Ah- C h) TOM ^ ' + ^ °tj ~ , 

(5.53) 

„ / 2roM(^0+2co) _ ( tz> i Tyf \ i 2thm(^H+2c h ) _ / r? i D/ ^ 

^ ~ V T -OM(Ao+2 C o)+THM(AH+2c H ) rHM ^ ffl + ^ Hi) + TOM (Ao+2 C o)+THM(AH+2 CH ) rOM ^ 0i + B Oi; 

(5.54) 

where r/j and 774 correspond to longitudinal and transverse acoustic phonons, respectively. 

6 The Kauzmann paradox and the VTF law; specific 
heat, relaxation times and transport coefficients 

First, we investigate the entropy due to intra-band density fluctuations. Since sound 
depends on the temperature, we must start with the interaction Hamiltonian. Here we 
consider the thermodynamical function. We introduce the Hamiltonian with a parameter 

A 0H3EEE1: 

H{\) = H + XHj. (6.1) 
The thermodynamical function Q\ for the above Hamiltonian is given by 

Vl x = ~ ln { Tre -/3(H () -/,iV+AH 7 )}_ ( 6 2 ) 
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Here we calculate the thermodynamical function due to intra-band density fluctuations 
denoted by flmtrax- From Eq.(2.13), we obtain 

= a >A= 2 ^ W/309 

= -^E^(g)^ A a (g), (6.3) 

where < • • • >\ and F x (q) are calculated by Hj replaced by XHj in < • • • > c and F(q). 
Replacing V a/ 3(q) by XVa^q) into Eq.(6.3), neglecting the A-dependence of P a (q) because 
of the condition f3h\ J M | <C 1 in Eq.(4.3), we obtain 

]TR 7 - Pa(g)A^(g)}F^(g) = P a (g)<f a/9 . (6.4) 



Eq.(6.4) leads to 



^oo ^ohA 1 /(l - XPhVhh)Pq XPqVorPh 
^ho f hh) det I |V APh^ho^h (1 - AP Voo)P H 



(6.5) 



det | | = (1 - AP Voo)(l - APhVhh) - A 2 P VohPhVW (6.6) 
Thus Eq.(6.3) yields 

^Wa 1^ (1- AP h Vhh)PqK)o + (1 - AP Vbo)PH Vhh + 2XP P H V 0U V U0 
d\ 2(3^ (1 - AP Vbo)(l - APhV^h) - A 2 P P h K>h^ho 

= ^J2§x ln{(1 " AP ° y ° o)(1 " APhV ™) " A'PoPhVohVho} (6.7) 

i 

Integrating Eq.(6.7) in A and putting A — > 1, we obtain the thermodynamical function 
due to the intra-band density fluctuations, VL nitvai : 

fiintra = ^ ^ M(l " WooXl " P^lffl) " PoPh^OH^Ho} + ^intraO, (6.8) 

where the term f2i ntr ao comes from the initial condition, In order to calculate the entropy 
we neglect l/r terms under the condition 13%/tq <C 1 and we put q — > \v n . From 
Eq.(4.3), P a (q) = (5{^ q f /(iu 2 n - «) 2 ), we obtain 

J_^, (^-^o(g))(^-^H(g))-/3 2 W 2 K^ 2 VbHVHO 

mtra 2^ 111 (i^-K ? ) 2 )K- W) mtra0 

_J_v-i (^ 2 -^ 2 +(g))(^ 2 -^s 2 -(g)) n , fiQ x 
" 2^ ^ 111 _ ( ^) 2 )K - K,) 2 ) intra0 ( j 
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The Contour integration of Eq.(6.9) gives 



fiintra = \ - e~^^){l - e^-M) - ln(l - e -«,)(l " e"^)] + ^intraO- 

(6.10) 

^intraO = ^ln(l-e-^<), (6.11) 
aqr 

where we have neglected the term | Ylqi^sq — ^°)- The first, the second and the third 
terms in Eq.(6.10) correspond to the thermodynamical function for the sound, the fluc- 
tuations due to the bubbles of the intra-band elementary excitations and the dissipation 
due to the diffusion, respectively. The third term compensates the second term and the 
system becomes a local equilibrium. From the instability of sound, Eq.(4.27), the dom- 
inant contribution of the temperature dependence to the thermodynamical function in 
the first and the second terms comes from the regions q ~ K. We separate the region 
of wavevector q into the number of states N near K and the remaining part. Using 
f3huj s ±(q), I3huj° aq <C 1 under the condition f3h\ J M | <C 1, we obtain 

i4 intra — p \ / T " sound 




~r ^ ' sound 



0M(1 + PVoo(K)(l +(3V HH (K)) -(3 2 V OK (K)V KO (K)} + n' sound , (6.12) 

fi'sound = ^£lHl ~ e-^W) + ln(l - e -^-W)], (6.13) 

where the prime on fi sound and £ means that the region of the wavevector q is limited 
in the remaining region with the number of states N — No excluded the region with 
the number of states N near K. Now we investigate the first term of Eq.(6.12) using 
Eq.(4.27). 

(1 + /3Voo(K))(l + pVjm{.K)) - t3 2 V K{K)V no {K) 

= {1 _ + To-iVooW + Vjmik) 

^(l-i)(2- ^ ) ;^ ) ) (6,4) 

Thus we obtain 

fiintra = ^ hl{l - ^} + fi'sound, (6.15) 

where we have neglected the constant term, which does not contribute to entropy. Note 
that the dissipation compensates the fluctuation entropy in the remaining region with the 
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number of states N — N . The entropy due to the intra-band density fluctuations, Si n t ra 
is given by 

Sintra = S K + S' sound , (6.16) 

^ _ NokB_r^ (617) 

2 T-Tn K ' 



s m , * - fcB yjV„(i-e-^ + <«)) + i„(i- e -^-(n + iE'{ 55 S^+ ii S^}- 

(6.18) 

Sk manifests the Kauzmann crisis. We call Sk the Kauzmann entropy. 

A state of iV-particles distributed randomly in space corresponds to a minimum of N- 
particle potential in a configuration space. A hopping of a particle from a site to a vacancy 
corresponds to a jump from a deep valley to another deep valley in the multi-dimentional 
configuration space in the energy landscape model (ELM) as the a-relaxation process. 
The successive hoppings constitute a configuration space. A hopping also generates intra- 
band elementary fluctuations and the successive hoppings yield the fluctuation entropy 
Sk due to the intra-band density fluctuations. The hopping probability is proportional 
to the configuration number, which is e SK ^ kB from the Einstein relation. The hopping 
probability of a particle is given uniquely by the hopping amplitude: 



zS 



K_ 



J = e Nk B = e T - T o, (6.19) 

where E = zN T /2N and z is of the order of the number of the surrounding H 2 
molecules. This equation manifests the Vogel-Tamman-Fulcher (VTF) law. The sound 
ou s ± consists of the individual sound frequencies u>o(q) and oun(q). As seen in the first 
terms on the right hand side of the first line in Eq.(6.12), the mixing between the entropy 
of sound u s ± and the fluctuation entropy due to the intra-band elementary excitations o>o q 
and u; Hq yields a unique hopping amplitude J. This fact originates from the temperature- 
dependent sound. Thus J governs all individual quantities v a and J aix . Then we can 
put 

\J afl \ oc J, (6.20) 



'U°T oc J. (6.21) 

Since all J afM are governed by a unique J, we can regard z as the order of the number 
of the surrounding H 2 molecules. It should be noted that the hopping J generates the 
intra-band elementary excitations, the sound u s ±(q) and the Kauzmann entropy Sk, while 
Sk determines J in Eq.(6.19). Thus u 8 ±(q) in Eq.(4.23) is self-consistently determined. 

Next we investigate the entropy due to inter-band density fluctuations. The thermo- 
dynamical function due to the inter-band density fluctuations is given in a similar manner 
to the sound: 

o J_y^ K-<(g))K-^-(<?)) , 
inter "2^ ln (K-^ )(K-^) +i2intcr0 ' 

= I ^[ln(l - e -^A+(«))( 1 _ e -/M*, A _(9)) _ ln(l - e -^°){\ - e-^ H )] + fi intcr0 , 

(6.22) 
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^intcro = \Y, Xn{l ~ e^X 1 - e ~^ H ) (6-23) 

' Aq 

where we have neglected the term | Y2\ q (hw\q _ ^) • The entropy due to the inter-band 
density fluctuations is given by 

Winter Sphonon 

+ S a - S n , S a = Sn + S n „, (6.24) 
W = -fa £[ln(l- e -^>) + ln(l-e-^-<*>)] + i ^[-ig^+ ^W ^ 

Aq Aq 

(6.25) 

Aq Aq 

In Eq.(6.24), the first, the second and the third terms correspond to the entropies for 
the phonons, the fluctuations due to the bubbles of inter-band elementary fluctuations 
and the dissipation due to the viscosity, respectively. In this case compared with the 
Kauzmann entropy, the phonon entropy and the fluctuation entropy do not mix. This fact 
originates from the temperature-independent phonons. The fluctuation entropy lowers the 
equilibrium entropy, but the dissipative entropy compensates the fluctuation entropy and 
the system becomes to a local equilibrium. 

A propagation of an up and down transition at a site to a surrounding site corre- 
sponds to a jump of a shallow valley to another shallow valley in the multi-dimentional 
configuration space in the ELM as the /3-relaxation process. The successive propagations 
constitute another configuration space different from hoppings. The probability of the 
magnitude of the randomness of harmonic frequencies of a-particles is proportional to the 
configuration number e Sfla ^ kB . The probability of the magnitude of randomness of the 
harmonic frequency of an a-particle is proportional to 

na = e * a SnJN kB = exp{-^^}, (6.27) 

where z a is of the order of the a- surrounding particles of an a-particles. Thus we obtain 

V / CC" 7 ocQ Q . (6.28) 

The Kauzmann entropy S'k diverges negatively at T so that the system seems to occur 
the entropy crisis. But the inter-band fluctuation entropy Sq crosses S'k above T and 
prevents the crisis. The Kauzmann entropy Sk dominates above the crossover temper- 
ature, but the fluctuation entropy Sn dominates below the crossover temperature. But 
below the crossover temperature, the dissipative entropy compensates the Sq completely. 
We can identify the crossover temperature with the liquid-glass transition temperature 



T ■ 

± 9- 



S K = S n \ T=T . (6.29) 



Sound is a collision wave essential in a fluid, while phonons are elastic waves essential in 
a solid. The glass transion is a sort of dynamical transition. 
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Next we investigate the specific heat. First, we calculate the specific heat due to the 
entropy of the intra-band density fluctuations, which consists of the Kauzmann entropy 
Sk and the entropy due to the conventional sound, S^und- The specific heat due to Sk, 
Ck is given by 

f Npk B T Tp r rp rp 

Ck = \ 2 {T - To)2 V m rt ■ (6-30) 
\ for T<T g K ' 

The specific heat due to the sound, C' sound is given by 

C - d ~ T\ [ si nh 2 ^±M + sinh 2 ^ = kB(N ~ No) ' (6 ' 31) 

where we have used the relation f3hui s ±(q) <C 1. C' sound remains ks(N — N ): constant 
around the liquid-glass transitionT s . Eqs. (6.30,31) shows the gap of specific heat at T g . 

Next we investigate the specific heat due to the inter-band density fluctuations. In 
this case, since the inter-band fluctuation entropy cancels with the dissipative entropy, 
there remains only the specific heat of phonons. This is because the structure of the 
phonons does change little above and below the glass transition, while the structure of 
the sound essentially depends on the temperature. We obtain the conventional specific 
heat of phonons, C p h non as 

Ph0n ° n 4 ^ sinh2 ^tM + sinh2 J - ^ 6Z) 

At low temperatures, C p h non oc T 3 . 

Now, we investigate the VTF law on the transport coefficients, the relaxation times 
and the velocities of the modes. 

(i) T < T: 

The velocities of the particles and the sound only depend on J a0 , so we obtain 

v ap , v s ± oc J, v\± = constant. (6.33) 

Here it should be noted the phonon velocities v\± is constant. 

(ii) T <T g < T: 

The term \/U a ^ v oc J afl dominates; the a-relaxation. We obtain 

T a0, T aM « J (6.34) 

If we consider Eqs. (6.33), (6.34) and diffusivity in (4.30) and viscosity r] X in Eqs. (5.53) 
and (5.54), we obtain 

D ± , rj^ 1 oc J. (6.35) 

It should be noted that the sound velocity v s ± oc J in Eq.(6.33) plays an essential role in 
the diffusivity D ± . Eq.(6.35) satisfies the Stokes law. 

(iii) T < T < T g : 

The term s/TTTT™ oc Vt a dominates; the /3-relaxation. Thus we obtain 

T «0> T aM « H, (6.36) 

D ± oc J 2 , r] X oc SI' 1 . (6.37) 
In this regime there remains little diffusion. 
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7 Concluding Remarks 



We have calculated the sound and the diffusion from intra-band density fluctuations, 
and the phonons and the viscosity from inter-band density fluctuations in water taking 
into account of the interaction potentials Vbo, Vhh and Voh, where we have included the 
pair distribution functions goo, gun, and gon- The terms Voh consists of Vq^ and Vq H , 
where the pair distribution function of which are given by g Q c ^ and g Q \, respectively. 
In water the hydrogen bonding plays an essential role in constructiong the hydrogen 
bonding network. But in the nearest neighbour approximation the < VohVho >c includes 
the unphysical terms such as Vq^V^q' and Vq^V^q. Since we concentrate ourselves to 
the hydrogen bonding network, we have taken into account of only the physical terms, 
^OH b ^HO an d ^oh^ho 5 ' navm g comfirmed that when the term VohVho appears, VohVho 
should be replaced by V^V^ + Vg&V^. 

In water the interaction potential Vqh and the pair distribution function g^ are the 
most important. The peak of g \ relates to the hydrogen bonding, which clusters water 
molecules. Thus the sound instability occurs at a reciprocal particle distance correspond- 
ing to a hydrogen bond length K = Kq\. If the magnitude of Vo H (a b H ) is large enough 
compared with the magnitude Vbo(aoo) and Vhh(ohh), the uj s -(q) and v s - in Eqs.(4.23), 
(4.24), and the diffusion coefficient D~ in Eq.(4.30) disappear. 

In the term Vhh we have taken into account of the interaction between hydrogens in 
different water molecules so that the term relates to the hopping of hydrogens. Since 
we neglect the interaction between hydrogens in the same water molecule counting the 
number of hydrogens, 2N, in phonons we obtain the 3-acoustic modes, and the 3-optical 
modes which are doubly degenerate. Since the number of hydrogens in a water molecules 
is 2, if we include the interaction between hydrogens in the same water molecule, the 
degenerate 3-optical phonon modes split into 6-optical ones. 

We have elucidated the Kauzmann paradox on the entropy crisis which originates 
from the sound instability at a reciprocal particle distance K = Kq\ corresponding to a 
hydrogen bond length and at an instability temperature T . The Kauzmann entropy Sk 
originates from the mixing between the sound and the intra-band fluctuation entropies 
which determine a unique hopping amplitude J. J governs the velocities v a and the 
hopping amplitudes J^ a for the individual particles and sound velocity. But the relative 
magnitude of the individual particle velocities are determined by the coupling strength in 
Eq.(2.8). On the other hand, since the mixing between the phonon and inter-band fluc- 
tuation entropies does not occur, the inter-band fluctuation entropies for the individual 
particles Sn a are independent. Thus the inter-band fluctuation entropy Sq = Sq q + Sq h 
prevents the Kauzmann entropy crisis at T g . Below T g the remaining intra-band fluctua- 
tion entropy free from the mixing and the inter-band fluctuation entropy are compensated 
by the dissipative entropy due to diffusion and viscosity, respectively; there remain the 
sound entropy only with the remaining number of states N — N free from the mixing, 
and the complete phonon entropy. 

We have elucidated the VTF law on the relaxation times and the transport coefficients 
taking into account of the random scattering processes due to the random hopping ampli- 
tudes and the random harmonic frequencies. The random hopping amplitudes determine 
the VTF law and the a-relaxation, while the random harmonic frequencies determine the 
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/^-relaxation. 

We have calculated phonons and viscosity taking into account of the electric polar- 
ization. Since the electric poralization relates to the p-states of particles, the electric 
polarization connects the inter-band density fluctuations and does not connects the intra- 
band density fluctuations. The effect of the electric polarization on phonon modes is to 
separate the longitudinal and transverse modes constant values according to the param- 
eters c a in Eq.(5.41). We have also elucidated the boson peaks Eqs.(5.49, 50). 
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